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ABSTRACT 



Context. The light elements Li, O, Na, Al, and Mg are known to show star to star variations in the globular clusters 47 Tuc and 
NGC 6752. Such variations are interpreted as due to processing in a previous generation of stars. 

Aims. In this paper we investigate the abundances of the a-element sulphur, for which no previous measurements exist. In fact this 
element has not been investigated in any Galactic globular cluster so far. The only globular cluster for which such measurements 
are available is Terzan 7, which belongs to the Sgr dSph. 

Methods. We use high resolution spectra of the Si Mult. I, acquired with the UVES spectrograph at the 8.2m VLT-Kueyen 
telescope, for turn-off and giant stars in the two globular clusters. The spectra are analysed making use of ATLAS static plane 
parallel model atmospheres and SYNTHE spectrum synthesis. We also compute 3D corrections from CO^BOLD hydrodynamic 
models and apply corrections due to NLTE effects taken from the literature. 

Results. In the cluster NGC 6752 sulphur has been measured only in four subgiant stars. We find no significant star to star scatter 
and a mean {[S/Fe]> = +0.49 ± 0.15, consistent with what observed in field stars of the same metallicity. In the cluster 47 Tuc we 
measured S in 4 turn-off and 5 subgiant stars with a mean ([S/Fe]) = +0.18 + 0.14. While this result is compatible with no star 
to star scatter we notice a statistically significant correlation of the sulphur abundance with the sodium abundance and a tentative 
correlation with the silicon abundance. 

Conclusions. The sulphur ~ sodium correlation is not easily explained in terms of nucleosynthesis. An origin due to atomic 
diffusion can be easily dismissed. The correlation cannot be easily dismissed either, in view of its statistical significance, until 
better data for a larger number of stars is available. 
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1. Introduction 



Galaxy: globular clusters: individual: NGC 104 - Galaxy: globular clusters: individual: NGC 



' There is substantial consensus that the so-called a-elements (O, 
Ne, Mg, Si, S, Ar, Ca, and Ti) are synthesized essentially dur- 
ing hydrostatic and explosive burning phases in massive stars, 
and then released in the interstellar mediu m when such stars 
. undergo type II Supernova (SN) explosion jChieffi & Limongil 
■ 120041) . As such, a-elements provide crucial information to re- 
' construct the star formation history of stellar populations. O, 
Mg, Si, Ca, and Ti are easily measured in stellar spectra and are 
thus the most frequent choice in stellar abundance studies. S on 
the other hand is the natural choice when abundances are de- 
termined in the Interstellar Medium (ISM), Planetary Nebulae 
(PN), or extragalactic absorption systems (e.g. Damped Lyman 
Alpha systems, DLA) due to the ease of measuring it and its 
weak tendency to form dust. The extent to which stellar and 
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ISM / PN / DLA a-element enrichments are comparable is thus 
dependent on wether or not S is supposed to vary in lockstep 
with the other (T-elements. 

Sulphur abundances in Milky Way (MW) st ars have been 
obtained by several author s in th e latest years ( iNissen et alJ 
2007; Caffau et al. 20 071; iRvdd [2006; Caffau et al.ll2005b : 
Korn & Rvddi2005l; INissen et al. 2004; Rvde & Lambertll2004 . 
among the most recent), but agreement on its behavior at low 
metalliciti es is not reached y et. While some of the cited au- 
thors (e.g. INissen et al.ll2007l) find S to follow the same trend 
as other c-elements ([S/Fe] showing a plateau at low metallic- 
ities), others point to wards a steady growth of [ S/Fe] with de- 
creasing metallicity dlsraelian & Reboldl200lL Takada-Hidai 
et al. 20 021), or an increase in the scatter of [S/Fe] at low 
[Fe/H] dCaffau et al.ll2005bl) . The cited debate can in fact be 
extended to the overall behavior of a-elements; while the stan- 
dard view of [(T/Fe] evolution with metallicity has for some 
ti me now assumed th e presence of a plateau below [Fe/H] — 
1 (lMcWilliamlll"997l) . observations as well as models suggest 
the possibility of a continuing rise down to the lowest observed 



metallicit ies. at least for Mg and O (first pointed out bv A bia 
& Rebolo 1989, an d more recent ly bv lFrancois et al.ll2004L and 
lGonzale"z Hernandez et al.ll2008l) . 

It is perhaps surprising that sulphur has been studied up 
to now almost exclusively in field stars. To our knowledge, the 
only analysis of the S abundan ce in globular cluste rs is also the 
only one for extragalactic stars (ICaffau et al.l2005aL in Terzan 7, 
one of the globular clusters of the Sagittarius dwarf Spheroidal 
Galaxy), where it behaves consistently with the rest of the a- 
elements studied in the cluster, showing a low ratio with iron 
typical of dwarf sphe roidal systems ( Tautvai siene et al.ll2004l : 
ISbordone et al.ll2007h . But on the other hand, S measurements 
are totally missing in MW Globular Clusters (GC). Globular 
Clusters are somewhat peculiar environments regarding chemi- 
cal enrichment, displaying some peculiar abundance anomalies 
and abundance correlati ons which have neve r been found else- 
where (for a review see lGratton et al.ll2004l) . Among them the 
most well known is likely the so-called O-Na anticorrelation, 
which appears to be widespread in MW globular cluster stars, 
and the Mg-Al anticorrelation. The nucleosynthetic processes 
responsible for this are believed to be essentially understood, 
since the chemical yields match well the ones derived from p- 
capture processes occurring in (or near) regions where H is be- 
ing burned via the C NO cycle at sufficiently high temperatures 
(iGratton et al.H2004l) . On the other hand, such a condition is 
reached during different phases of stellar evolution, so that the 
actual processing site is still debated, as well as the mechanism 
through which the enriched gas is dispersed inside the clus- 
ter (^n^ura_e^_al^ 2002; D'Antona ej^al. 2005; Decressin et 
al. l2007h . There is nevertheless some consensus that the yields 
should be released through slow moving ejecta, thus allowing 
the GC gravitational potential to retain them more easily. 

We present here the first sulphur abundances in the GC 47 
Tuc and NGC 6752, both of which display evident Na-O and 
Mg-Al abundance anticorrelations. They are also the only two 
clusters where a Na-Li anticorre lation has also been observed 
(IPasquini et"ani2005l;lBonifacio et al., 2007.) . 

2. Observational Data 

The spectra analysed in this paper have been acquired 
in the course of the ESO-Large Programme 165.L-0263. 
Observations, atmospheric parameters an d abundances for ele - 
ments other tha n S hav e been presente d inlGratton et al.l (1200 ih . 
ICarretta etal] ( l2004l) . and Bonifaci o~et alj ( 120071) We are 
here concerne d only with th e dichroic #2, CD4 cross-disperser 
UVES@VLT dPekker et al.i.200 0) spectra (setting centered at 
760 nm). For most observations the slit width was of 1" for a 
resolution of R~43000. For the purpose of the sulphur mea- 
surement, only upper-red CCD data were used. Typical Signal- 
to-Noise ratios (S/N) in the range were of about 15-20 per 
resolution element. Sulphur abundances were determined for 
nine stars observed in 47 Tuc, four of them being turn off' (TO) 
stars, 5 subgiants (SG). Four stars, all subgiants, were con- 
sidered in NGC 6752. Typical exposure times were Ih per star 
for NGC 6752 SG stars, 2h per stai's for 47 Tuc SG and 4h per 
stars (in both cases performing multiple Ih exposures) for 47 
Tuc TO stars. An example of S i mult. 1 region is shown in fig. 

m 




921.0 921.5 922.0 922.5 923.0 923.5 924.0 924.5 
wavelength (nm) 

Fig. 1. The (normalized) spectrum of 47 Tuc star # 1012 
with superimposed a synthetic spectrum with T^ff 5832 K, 
log^ 4.05, [Fe/H]=-0.5. The three fines of Si mult. 1 are 
labeled, as well as the most prominent telluric absorptions 
(marked with © symbol). 



3. Analysis 

The atmospheric parameters have been taken fr om Gratton et 
al. (2001) for the stai's in NGC 6 752 and from Carretta et alJ 
(2004]) and lBonifacio et alj ( l2007l) (for star # 952) for the stars 
in 47 Tuc. 47 Tuc stars have a typical temperature of 5800 
K for the TO stars (log g=4.05) and of 5100 K for the sub- 
giants (log g=3.84), the four NGC 6752 share the same atmo- 
sphere parameters, T^s -5347 K log g=3.54. For each star we 
comp uted a rnodel atmosphe re using version 9 of t he ATLAS 
code (Kurucz" 1993b', '2005a) running under Linux ("Sbor donj 
2005; Sbordone et al. 2004 ). We used the up dated Opacity 
Distribution Functions o f Castelli & Kuruc3 (12003 ) with micro- 
turbulent velocity of 1 kms ' and enhancement of a-elements. 
The synthetic spectra w ere computed using th e SYNTHE suite 
("Kurucz" 1 993b, '2005a|) running under Linux (ISbordonell20"05l ; 



^bordone et al. 2004 

The strongest lines of S i in the optical range are the lines 
of Mult. 1 (three lines around 923 nm) and of Mult. 6 (two lines 
at 869 nm, see Table[T]). We were unable to convincingly detect 
the lines of Mult. 6 in any of our spectra, we therefore did not 
even search for Mult. 8 ( around 675 nm) which is weaker. Our 
analysis is therefore based exclusively on the lines of Mult. 1 . 
We adopt the log^/ value s of the NIST database, which holds 
the values of IWiese et al.l d 19691) . which are experimental but 
of D quality, corresponding to a possible error up to 50 % (see 
Table[T]for the values). A thorough discussion of the ava ilable 
values for S i atomic data is avilable lCaffau et aTl ( l2005bl) . 

The lines are in a spectral range contaminated by telluric 
absorption, so that they are rarely all free from telluric contam- 
ination. Moreover, due to the low metallicity of the two clus- 
ters examined and to the low S/N ratio of the spectra, in most 
cases, only one of the S i lines of Mult. 1 is detectable. Due 
to the radial velocity of the clusters, the line which happens to 



Table 1. Sulphur line-by-line abundances for the program stars. 
The wavelengths and log gf if the three lines are indicated as 
well. 



A 


921.286 nm 


922.809 nm 


923.754 nm 


log gf 


0.42 


0.26 


0.04 


NGC 6752 


1406 


5.80 


6.07 


6.04 


1665 




6.35 


6.29 


1461 




6.35 


6.37 


1481 




6.35 


6.36 


47 Tuc 


1012 




6.56 




1081 




6.68 




952 




6.63 




975 




6.72 




201075 


6.65 


6.86 


6.43 


206415 


6.73 






429 




6.87 




433 




6.87 




482 




6.80 





be the most often unblended is the one at 922.8 nm, that is the 
line laying on the wing of the hydrogen Paschen f line. The 
measurement of the equivalent width (EW) of lines falling on 
the wings of broad features is not easy, and produces a value 
whi ch is usually lower , with respect to the true EW, by about 
5% (ICaffau et al.l l2005b1. To determine the sulphur abundance 
a line profile fitting procedure, taking into account the profile 
of the Paschen ^ is to be preferred. The other two lines of Mult. 
1 are far enough from the Paschen ^ not to be affected by it, 
so that for these lines the EW measurement would be a good 
choice. To keep the analysis as homogeneous as possible we 
nevertheless decided to determine the sulphur abundance us- 
ing line profile fitting for all the unblended features. Such fit- 
ting was performed bv means o f the code described in Caffau 
et al. (l2005bfK In the fitting procedure, the Une broadening 
was kept fixed at the instrumental level (7 km/s). At this res- 
olution, and with this S/N ratio, typical macroturbulent and/or 
rotational broadening would not affect the measurement sig- 
nificantly. We discarded the lines contaminated by telluric ab- 
sorption because the telluric line subtraction from a spectrum 
decreases the S/N that is already low. Also, we did not pos- 
sess fast rotator spectra taken at the same time as the science 
exposures, which would have imposed to rescale the telluric 
features, an operation which would add further uncertainty to 
the measurement. Line-by-line S abundances are listed in Table 

m 

The S I lines of Mult. 1 are known to suffer from departures 
from LTE jTakeda et alj|2005l) . We the refore appUed the NLTE 
corrections, interpolated in table 2 of iTakeda et alj (1200 5) to 
take this effect into account. For our programme stars these are 
all negative and range between -0.1 and -0.3 dex. 

3.1. Estimation of analysis uncertainties 

Uncertainty in the S abundance determination stems from both 
the limited S/N ratio of the spectra and the systematics of the 
analysis. To determine the effect of S/N we followed the same 



method described in ICaff^au et al.l (l2005bl) . namely by perform- 
ing Monte Carlo simulations. In that paper the error estimates 
are obtained for a sample of dwarf stars with [Fe/H]=-1.5. We 
have performed additional Monte Carlo tests to assess the Teff 
sensitivity of the S lines in the SG stars (rgft being the only pa- 
rameter showing a significant difference between TO and SG 
stars), and to ascertain that such uncertainties apply to the more 
metal rich objects of 47 Tuc and c ame up with v alues which are 
essentially identical to the ones in ICaffauet al.l (l2005b), which 
we will then assume from now on. 

Concerning random errors, the result of a 1000-events sim- 
ulations, for a parameter set compatible with the TO stars set, 
indicates an Icr uncertainty (for a single line) of 0.12 dex for 
S/N=30, and 0.22 dex for S/N=15. The same values hold for 
the SG stars. Given the typical S/N of our spectra, we will as- 
sume henceforth a typical uncertainty of 0.2 dex to be due to 
the S/N. 

For the TO sample, a change in log g of +0.25 dex causes 
a + 0.08 dex variation in A(S), + 100 K in temperature cause 
a variation of +0.06 dex in A(S), while [Fe/H]+0.2 dex in- 
duces a ± 0.03 dex variation in A(S). Comparing with table 3 of 
ICarretta et al.l (|2004) one notices that both T^fi and log g varia- 
tions act in opposite senses on [Fe/H] and on A(S). On the one 
hand, log g has a maximum stated error in lCarretta et al.l (l2004l) 
of 0.08 dex. Taking into account both the effect on metalUcity 
and the one on A(S), this would lead to a variation of 0.03 dex 
in [S/Fe], negligible in the present scope. 

On the other hand, ICarretta et al.l (|2004|) claim a maximum 
Tetf error of + 90 K for TO stars and +60 K for subgiants. 
Shifting the whole temperature scale rigidly would not make 
or break the [S/F e] spread we observe, but one could argue that 
the ICarretta et a l. ( 2004) Teff scale might be excessively "com- 
pressed". By lowering the TO stars temperature and rising the 
one of the subgiants by the maximum amount stated, one would 
reduce significantly the observed [S/Fe] spread. A 90 K lower 
temperature for the TO stars would lower their [Fe/H] by 0.08 
dex. The effect on S would then be due both to the reduced Tgff 
and to the decreased metallicity, the total effect being roughly 
AA(S)=+0.04 dex. A similar calculation for the SG stars leads 
to AA(S)=-0.03, which become A[S/Fe]=+0. 12 for the TO and 
-0.09 for the SG. This would reduce the [S/Fe] spread by 0.21 
dex, eliminating something less than half of it. On the other 
hand, such a change would lead [Fe/H] to differ between TO 
and SG stars. Currently, the two subgroups give precisely the 
same metallicity (-0.67) with a dispersion of 0.02 dex among 
TO stars and 0.01 dex among SG stars. Applying the afore- 
mentioned opposite temperature differences would separate the 
TO and SG stars by 0.14 dex in metallicity, roughly five times 
the quadratic sum of the two group uncertainties. We there- 
fore suggest that any Teff bias is very likely to be a simple 
additive shift to both TO and SG stars temperature. It is also 
worth noting that some imbalance (~0.13 dex) exists between 
Fe I an d Fe n abundances in the giants sample of ICarretta et al.l 
(120041) . According to Carretta and collaborators' error budget, 
this would be entirely solved by rigidly increasing the Tgff in 
subgiants by 100 K. Again, a rigid shift of the Teff scale would 
not affect the spread in A(S) nor in [S/Fe], and would also break 
the ionization equilibrium in the TO sample. Shifting the giant 
sample only would reduce the spread in sulphur abundance, but 



Table 2. 3D corrections for the S i 922.8 nm line, computed 
for a set of stellar parameters compatible with the ones in our 
sample, for an EW of 12 pm. The corrections are in the form 
Abu(3D)-Abu(lD). The model parameters are listed as well. A 
microturbulence of 1.5 km s ' has been employed in the LHD 
ID model. 







[Fe/H] 


Corr. 


5000 


4.0 


0.0 


0.141 


5000 


4.0 


-1.0 


0.140 


4920 


3.5 


0.0 


0.146 


4930 


3.5 


-1.0 


0.167 



Table 3. Sulphur abundances for the s ubgiant stars in NGC 
6752. The star identification comes from Gratt on et al ] (l2001h . 
as well as the atmosphere parameters: T^ff =5347 K log g=3.54 
Vturb= 1-1 kms"' , [F e/H1=-1.43 for all sta rs. The abundances 
of Na and Mg are from iOratton et al.l (1200 ih . The number of S i 
lines used is given in the second column. 



Star 


N 


[S/Fe] 


[Na/Fe] 


[Mg/Fe] 


1406 


3 


+0.19 


+0.02 


+0.20 


1665 


2 


+0.54 


+0.10 


+0.10 


1461 


2 


+0.58 


+0.29 


+0.10 


1481 


2 


+0.57 


+0.54 


-0.07 



generate a Fe abundance spread. Besides, any increase in the 
subgiants Jeff would mainly affect F ei, raising its abund ance. 
A discrepancy already exists between Carretta et al. ( 2004) and 
[ Koch & McWilliam (2008) [Fe/H] scales, Koch & McWilliam 
(i2008h deriving an average [Fe/H] lower by 0.09 dex. Such a 
discrepancy would b e worsened, should the temperature for the 
ICarretta et alJ (l2004l) subgiants be increased. 

3.2. 3D corrections 

We computed 3D corrections for the S i 922.8 nm line, for 
a set of atmosphere models encompassing the parameters of 
our program stars. In Table |2] we present the resulting values. 
Assuming an equivalent width of 12 pm (typical for the ob- 
served stars), we tabulate the difference between the abundance 
needed to reprodu ce the E quivalent Width (EW) in the 3D 
CO^BOLD model dFrevtaget al. 2002; Wedemeyer et al. 2004) 
and the one needed when using the comparison IDlhd model 
(ICaffau et al.ll2007h . As such, the 3D abundance is obtained by 
adding the correction to our ID abundances. Table |2] reports 
the model parameters as well. 3D corrections appear to be on 
average 0. 15 dex and positive, and appear highly homogeneous 
across the explored parameter space. In particular, the differ- 
ences in the corrections for TO and SG stars are small. They 
have not been added to our measurements in the table and in 
the figures, since comparison elements, in this and other stud- 
ies do not have 3D corrections applied. 

4. Results 

Our results are summarised in Tables [3] and |4] For the reader's 
information we report there also some information on our pro- 
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Fig. 2. [S/Fe] plotted against [Na/Fe] for our sample of 47 
Tuc stars. Na spread is added for reference from Brown & 
Wallerstein (119921) (open red triangles), Alves-Brito e t alJ 
(2005) (open green stars), and Koch & M cWilliam (2008h 
(open blue squares). A typical error for our dataset is also dis- 
played. Stars in our sample belonging to the TO are identified 
by a larger open circle around them. 



gr am stars extracted fromlOratton et al.l (1200 ih and Carretta et 
al. (|2004[) . The other abundances coming from these papers are 
always ID-LTE values, with the exception of N a, wh i ch ha d 
NLTE corrections applied according to Gratton et al.l (Il999l) . 
We adopted a solar sulphur abundance of A(S)=7.21 in the 
[S/Fe] determination. 

In NGC 6752, due to its low metallicity ([Fe/H]=-1.43), 
we were only able to measure the sulphur abundance in the 
four subgiant stars. In the fainter TO stars the S/N ratio was 
too low to detect the S i lines. From the four subgiants we de- 
rive < [S/Fe] >= +0.49 + 0.15. Our estimated error on the 
sulphur abundances is 0.2 dex, dominated by the noise in the 
observed spectra. Thus the dispersion around the mean is fully 
compatible with the errors. Moreover, [S/Fe] in NGC 6752 is 
fully compatibl e with the value fou nd among field stars at this 
metallicity (see lCaffau et al.ll2005b figure 10). 

For the more metal-rich cluster, 47 Tuc, we have been able 
to measure at least one S i line for all the stars. In star 201075 
all the three lines of Mult. 1 were measurable. The [S/Fe] ratio 
varies by almost 0.5 dex, with an average value of < [S/Fe] >= 
0.18 + 0.14. Given that the uncertainty in the S abundance in 
each star is estimated at 0.2 dex, the observed spread can in fact 
simply reflect the measurement uncertainty. Moreover, as in the 
case of NGC 6752, the mean [S/Fe] is totally compatible with 
the valu e displayed by field stars of abundance comparable to 
47 Tuc (ICaffau et al.ll2005bl) . 



Table 4. Sulphur abundances for the stars in 47 Tuc. The star numbers come from Carretta et al.l (|2004 ). as well well as atmo- 
spheric parameters and abundances of Na, Mg and Si. For star # 952 see lBonifacio et all (l2007l) . The number of S i lines used is 
given in the fourth column. 



Star 




logg N 


[S/Fe] 


[Fe/H] 


[Na/Fe] 


[Mg/Fe] 


[Si/Fe] 




K 


cgs 


dex 


dex 


dex 


dex 


dex 


1012 


5832 


4.05 1 


-0.01 


-0.64 


-0.14 


+0.45 


+0.03 


1081 


5832 


4.05 1 


+0.11 


-0.64 


-0.34 


+0.50 


+0.13 


952 


5832 


4.05 1 


+0.06 


-0.64 


-0.08 






975 


5832 


4.05 1 


+0.15 


-0.64 


+0.22 


+0.21 


-0.10 


201075 


5165 


3.84 3 


+0.12 


-0.68 


+0.11 


+0.42 


+0.38 


206415 


5112 


3.84 1 


+0.22 


-0.70 


+0.10 


+0.34 


+0.45 


429 


5081 


3.84 1 


+0.31 


-0.65 


+0.31 


+0.39 


+0.42 


433 


5106 


3.84 1 


+0.44 


-0.78 


+0.24 


+0.49 


+0.21 


482 


5090 


3.84 1 


+0.21 


-0.62 


+0.06 


+0.64 


+0.28 



0.6 



0.4 



0.2 



0.0 



-0.2 



^1081 

,g,1012 

A AAA 



-0.2 



0.0 



0.2 
[Si/Fe] 



0.4 



0.6 



Fig. 3. [S/Fe] plotted against [Si/Fe] for our sample of 47 
Tuc stars. Si spr ead is added for reference from Brown & 
Waller stein (il992h (open red t riangles), ves-Brito et alj 
(l2005h (open green stars), and iKoch & McWilHam (2008) 
(open blue squares). A typical error for our dataset is also dis- 
played. Stars in our sample belonging to the TO are identified 
by a larger open circle around them. 



47 Tuc is believed to display Na-O as well a s Mg-Al abun- 
dance anti-correlations dCarretta et alj 12004*). We thus de- 
cided to check for possible correlations between S and some 
of the cited elements. Sulphur has never been measured be- 
fore in 47 Tuc but interesting insight can be gained regarding 
other elements by comparing Carretta et al. (2004) with the 
ones obtained bv Brown & Wallerstein ( 1992), Alves-Brito et 
al.12005), and Ko ch & McWilUam (.2008 ). In Fig. |2] [S/Fe] is 
plotted against [Na/Fe]. Black filled circles represent our data, 
points surrounded by a slightly larger circle indicate TO stars. 



and every data point is flagged with the corresponding star 
number. The tvpical size for the error is also indicated. Brown 
& Waflerstein (1992), Alves-Brito et al. (2005), and Koch & 
Mc William (l2008l) do not measure sulphur, but for reference 
we add the ir values for [Na/Fe l: blue open squares for Ko ch & 
Mc Wi lliam (l2008h red open triangles for Brown & Wallers teinl 
(Il992h . and green open stars for Alves-Brit o et al.i (|2005). The 
correlation between Na and S abundance is clearly apparent in 
our data, and a Kendall t rank correlation test in fact deter- 
mines that the likelihood of [Na/Fe] and [S/Fe] to be correlated 
is 97. 8 "^B Another interestin g thing to notice is how all the 
studies with the exception of Koch & M cWilliamI (l2008l) detect 
a significant spread in the Na abundance. 

In Fig. [3] we instead plot [S/Fe] versus [Si/Fe], all the sym- 
bols being the same as in Fig.|2] Again, there is an evident sug- 
gestion of a correlation, with the exception of star 433 which, 
while having a very high S abundance, presents a less extreme 
[Si/Fe]. Likelihood of the correlation from Kendall r test is 
here about 79%, but grows to 91% if star 433 is removed from 
the sample. By removing star 433 from the [Na/Fe] vs [S/Fe] 
correlation test, one gets on the other hand a slight decrease 
of the correlation likelihood, down to 94%. It is worth noting 
how star 433 has, according to lCarretta et al.l (12004 ). an unusu- 
ally low [Fe/H] (0.14 dex below the average, way beyond the 
scatter of the Fe i abundances), and an unsatisfactory Fe i-Fe ii 
ionization equilibrium (0.21 dex difference). Since no evidence 
ever emerged of 47 Tuc showing a spread in iron abundance, 
one might regard these results as suggestive of some problem in 
the analysis of star 433. If we tentatively assign the average Fe 
abundance ([Fe/H]=-0.66) to star 433, we obtain [S/Fe]=0.32, 
more in line with the rest of the sample. In fact, star 433 has a 
[S/H] exceeding by just 0.14 dex the average value (([S/H])=- 
0.48. [S/H]433=-0.34). 

Given the fact that most low-S stars (and low-Na stars) 
are TO stars, it is worth investigating whether diffusion in 
the stellar atmosphere might be responsible for the observed 
abundance spread and its correlation with the Na abundance. 
One would in this case expect TO stars to be most affected 



' This is equivalent to say that the likelihood of the null hypothesis 
is 2.2 %. This is the likelihood that random errors applied to otherwise 
uncorrelated distributions might produce a level of correlation (tk) 
equal or superior to the observed one. 



by diffusion, while SG stars should not display them anymore 
due to the deeper reach of their atmospheric convective zone. 
According to current models this is not a viable explan ation of 
the observed spreads. According to iKorn et alJ (l2007h (fig. 1), 
by using a model with radiative acceleration and turbulent mix- 
ing, an old [Fe/H] =-2 TO star shows photospheric depletion of 
about 0.25 dex in Na, 0.20 dex in S and 0.15 dex in Fe. The 
effect would be then 0.05 dex on [S/Fe] and 0. 1 dex in [Na/Fe], 
too small to explain the observations. 

5. Discussion 

The sulphur abundance in NGC 6752 agrees with the general 
trend of a-elements in the field, and does not present any signif- 
icant spread. In 47 Tuc the situation is less clear. The presence 
of a significant correlation between the S and Na abundance, 
and the hint of a S/Si correlation, suggest that the S abundance 
spread might indeed be real. 

Recendy iKoch & McWilliamI (l2008l) . in a thorough anal- 
ysis of 47 Tuc, based on high quality spectra of bright giants, 
found only a very limited spread in abundances, and were un- 
able to confirm the presence of a significant spread in O, Na, 
Al and Mg. Su ch discrepancy with the results of Brown & 
Wallerstein CI 992). 'Carretta et al.l (|2004 . lAlves-Brito et al.1 
(I2005h . andteonifacio et al. (200^ has not found to date a sat- 
isfactory explanation. This is not the venue to reassess in detail 
the work of the cited authors. Suffice to say that, applying an 
error analysis similar to the one above sketched for sulphur, 
the reality of the spreads in 47 Tuc as well as their correlation 
appear difficult to dismiss. 

In FiglHwe plot [Si/ Fe] versus fMg/Fe] for our stars (abun- 
dances a re taken from ICarretta et alJ |2004|), the sample of 
[Brown & WallersteinI (Il992h and the one of lKoch & McWilliamI 
62008). Symbols are the same as in Fig.|2]and[3] The first thing 
that can be noted here is that, with the exce ptio n of star 4820 , 
our sample and the one of iKoch & McWilliamI (l2008h present 
in fact the same spread in [Mg/Fe], but a very different one 
in [Si/Fe]. A word of caution should be here given regarding 
silicon: it appears that Si shows a systematic discrepancy be- 
tween TO_and^^BJiiexti^^ stars (Mc William 
119971: ICohen et al.ll2004l: iBonifacio et all 2009). The reasons 
of the discrepancy are not known, NLTE effects being a pos- 
sibility. The sense and amount of the discrepancy is compati- 
ble with the one seen here, so there exists the possibility that 
the Si spread is not real. On the other hand, such an effect has 
never been reported on stars as metal rich as the ones in 47 Tuc. 
Another element for which the various authors do not reach a 
satisfactory a greement is Ca. Average TCa/Fel is 0.34 in Koch 
& Mc William (2008), 0.18 in Carretta et alj |2004j), but solar 
ac cording to Brown & Wallerstein (1992) and Alves-Brito et 
al. (12005b . 

Another important aspect is of course constituted by the ex- 
istence of significant correlations among the abundances. One 
might consider as worrying the fact that all the TO stars of 
the sample exhibit a low S abundance. On the other hand, star 
201075, which is a subgiant, and the one star where all the three 
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^ Star 482 is indeed suspect in a mor e general way, prese nting the 
largest Fei - Feii imbalance of the whole lCarretta et al] ( l2004l) sample. 



Fig. 4. [Si/Fe] plotted a gainst [M g/Fe ] for ou r sample of 47 
Tuc sta rs, plus iBrown & Waller steinI (1 1992!) (open red tri- 
angles), Alves-Brito et al. (2005) (green stars), and Koch & 
Mc William (2008) (open blue squares). Stars in our sample be- 
longing to the TO are identified by a larger open circle around 
them. A typical error for our sample is also displayed. 



lines of Mult. 1 were measured, shows low S as well. By exam- 
ining Table [Hone might notice that the line at 922.8 nm gives 
the highest abundance for star 201075. Since this line is the one 
most frequently measured in the other stars, one might object 
that, had it been the only one measured, the star would have 
assigned an higher abundance, compatible with the other sub- 
giants. It has although to be noted that in NGC 6752 both this 
line and the one at 923.7 nm are systematically measured, and 
no sign exist of the latter giving a lower abundance. As detailed 
in sect. 13.11 the quoted Teff uncertainty would in principle be 
sufficient to remove the spread between TO and subgiant stars, 
if TO stars temperature was lowered and SGB stars tempera- 
ture increased. On the other hand, this would make star 201075 
even more discrepant, and induce a Fe abundance spread be- 
tween TO and SGB stars. 

Furthermore, a Monte Carlo simulation of 1000 events has 
been performed to check the likelihood that a significant cor- 
relation might emerge in the [S/Fe] vs. [Na/Fe] plot when the 
quoted errors are superimposed to a constant [S/Fe] for all the 
stars. Unsurprisingly, the result (~2% likelihood) confirms the 
result of the Kendall r null-hypothesis probability. It is very 
unlikely that random errors might have produced the observed 
Na / S correlation. 

In short, we believe that strong hints of a S abundance in- 
homogeneity in 47 Tuc do exist. It is nevertheless clear that the 
currently available data do not allow, due to a combination of 
low statistics and significant uncertainties, to clarify the mat- 



ter. Further study would thus be in order, specifically targeting 
a large number of stars from TO to RGB, and focusing on the 
elements that show significant spreads. Still, being unable to 
justify the current observations by means of either random or 
systematic uncertainties, we feel compelled to look for a nu- 
cleosynthetic explanation. The Na-0 anticorrelation is usually 
explained by self pollution within the cluster by slow Na-rich 
winds, either from AGB stars or from massive rotating stars. In 
general, the stars of pristine cluster composition are expected 
to show low Na abundance, while a higher abundance is indica- 
tion of "polluted" stars. Should Si and S abundance spreads be 
confirmed, and should they both correlate with the Na spread 
(which is quite uncertain in the Si case), this would pose at least 
two problems. 

First, since Na-poor stars are Si-poor and S-poor as well, 
the aforementioned model would indicate them as the "unpol- 
luted" stars. 47 Tuc would then have formed with an unusual 
low abundance of Si and S, at variance e.g. with Mg. Low 
[c/Fe] ratios are indeed observed in dwarf galaxies, but typ- 
ically attribu ted to environmental effects (e.g. Lanfranchi & 
Matteucci ,20()7ir and no evidence exist that they might priv- 
ilege one explosive product (e.g. Si) over another (e.g. Mg). 

Second, as just noted, both Si and S are believed to 
be, almost exclusively, products of SN II explosive burnings. 
Elementary models forbid globular clusters to be chemically 
enriched by internal SN explosions, since the total released en- 
ergy is comparable to the cluster binding energy, and even one 
single SN II explosion should remove all the available gas from 
the cluster potential well, and stop the star formation. This is 
the reason why "slow wind" mechanisms have been invoked to 
produce e.g. the Na pollution. At the same time, the lack of any 
detectable spread in Fe abundance in 47 Tuc constitutes an ad- 
ditional problem, since one would expect some amount of Fe 
to be produced by the SN along with a-elements. On the other 
hand, such simplistic considerations provide likely too crude a 
model of the actual phenomenon, since they assume a spheri- 
cally symmetrical, central explosion which energy is efficiently 
transferred to the cluster gas. SN II explosion are in fact rarely 
spherically symmetrical dWang & Wheeled2008h . Instead, they 
tend to release a large part of the energy along more or less 
collimated jets which are not efficient in transferring the en- 
ergy to the cluster ISM. Moreover, SN which explode in a non- 
spherically symmetrical way tend to produce larger fallback, 
which would in turn limit their Fe yield, thus eliminating the 
"missing Fe spread" problem. It is also suggestive that detailed 
models of stellar superclusters foresee a destructive interfer- 
ence among shock w aves from diffe rent supernovae, within a 
given critical radius dWiinsch et al.i r2008). As a result the in- 
ner part of the cluster would not only retain the yields, but also 
experience an accelerated star formation. Whether this might 
be the case on the smaller scale of normal GC (or whether 
some/all the GC might form as superclusters) is anyway uncer- 
tain. Core collapse SN are a source of Na as well, and we see 
that the Na abundance correlates well with the S abundance. 
The applicability of the aforementioned model thus relies on 
the (yet unexplored) possibility that high Na-yield and high a- 
yield (or at least high S-yield) happen together 

An alternate way to produce S and Si would be through 
proton capture reactions. This is the same source invoked to 



explain the Na overabundance, which would be an obvious 
advantage. On the other hand, S should this way be produced 
through the ^'P(p, ')')^^S reaction. The phosphorus abundance 
in 47 Tuc is unknown, and in general, knowledge about P 
production channels and abundances in stars are sketchy at 
best (see Caffau et al. 2007, and references therein). However, 
P in the Sun is roughly 3 dex less abundant than S. If the 
same abundance ratio holds at lower metallicities, it would be 
impossible to produce any significant amount of S at expenses 
of P. Phosphorus should then be formed at expenses of Si, 
which should then anticorrelate with S, which is not the case. 
The next abundant enough seed is Mg (Mg — > Al Si — > P 
S), which indeed might show weak signs of anti-correlation 
with Si (Fig. ID, although they are below any statistical signif- 
icance. Given the very scarce knowledge about P abundances 
and production mechanisms, one might even suppose the 
"original" 47 Tuc P abundance to be heavily enhanced, thus 
providing enough seeds to form ~0.2 dex of S. P abundance 
has never been measured in 47 Tuc, but its abundance has 
been determined very recently i n Horizontal Branch (HB) stars 
in NGC 6397 and NGC6752 dHubrig et alj|2009l) . where it 
appears strongly enhanced, by more than 2 dex in NGC 6752, 
and almost 3 dex in NGC 6397. Chemical anomalies in HB 
stars are usually attribut ed to diffusion effects, but the current 
models (iMichaud et al.l 12008) do not seem to be account for 
more than ~0.7 dex of P enhancement in the photospheres 
of HB stars of similar temperature. We thus think that the 
possibility of an enhancement of P abundance in globular 
clusters cannot be ruled out at the moment. 
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